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The Fermi-liquid theory of superconductivity is applicable to a broad range of systems that are candidates for 
unconventional pairing, e.g. heavy fermion, organic and cuprate superconductors. Ginzburg-Landau theory 
provides a link between the thermodynamic properties of these superconductors and Fermi- liquid theory. The 
multiple superconducting phases of UPt3 illustrate the role that is played by the Ginzburg-Landau theory in 
interpreting these novel superconductors. Fundamental differences between unconventional and conventional 
anisotropic superconductors are illustrated by the unique effects that impurities have on the low-temperature 
transport properties of unconventional superconductors. For special classes of unconventional superconductors 
the low-temperature transport coefficients are universal, i.e. independent of the impurity concentration and 
scattering phase shift. The existence of a universal limit depends on the symmetry of the order parameter 
and is achieved at low temperatures ksT <C 7 -C Ao, where 7 is the bandwidth of the impurity induced 
Andreev bound states. In the case of UPt3 thermal conductivity measurements favor an E\ g or Ei u ground 
state. Measurements at ultra-low temperatures should distinguish different pairing states. 



1. Introduction 

Theoretical investigation of unconventional pair- 
ing began with the publication by Anderson and 
Morel [Q, "Generalized Bardeen-Cooper-Schrieffer 
States and the Proposed Low- Temperature Phase of 
3 He". They studied the physical consequences of 
BCS pairing with non-zero angular momentum, in- 
cluding superfluid phases with spontaneously broken 
time-reversal symmetry When superfluidity was dis- 
covered Q in 3 He it was immediately clear that this 
was not a conventional s-wave BCS superfluid be- 
cause there was more than one superfluid phase. Fur- 
ther evidence for unconventional pairing came from 
many experimental results and shortly after their dis- 
covery the three superfluid phases of 3 He were undis- 
putedly identified as p-wave spin-triplet superfluids. 

The search for unconventional superconductivity, 
the metallic analog of superfluidity in 3 He, was given 
a boost by the discoveries of superconductivity in 
the class of heavy-fermion metals [|| , particularly in 
the U-based compounds of UBei 3 , UPt 3 , URu 2 Si2, 
UM2AI3, and UPd2ALj. Unusual temperature de- 
pendences of the heat capacity, penetration depth, 
and sound absorption led to conjectures that these 
materials were unconventional superconductors Q. 
Much more experimental information is now avail- 



able, and there is consensus that some heavy-fermion 
superconductors (if not all of them) show unconven- 
tional pairing. Interest in unconventional supercon- 
ductivity expanded further with reports of several 
experiments on cuprate superconductors that sup- 
ported earlier predictions Q of "d-wave pairing" . 

A rigorous classification of superconductors by 
the angular momentum of the Cooper pairs (i.e. s- 
wave, p-wave, d-wave pairing, etc.) is not appro- 
priate in crystalline materials because angular mo- 
mentum is not a good quantum number. How- 
ever, the terms "d-wave pairing", etc. are often 
used interchangeably with "unconventional pairing" 
for states in which the pairing amplitude sponta- 
neously breaks one or more symmetries of the crys- 
talline phase, i.e. 1Z * A(1Z * p^) ^ A(p/), where 
7Z e G = Qspin x Q space, x T x Ugauge represents an 
operation of the full symmetry group Q other than a 
pure gauge transformation, and 

A Q /3(P/) ~ (a P/Q a_p //3 ) , (1) 
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is the equal-time pair amplitude, or order parameter, 
that describes BCS-type superconductors. Note that 
Pf is the momentum of a quasiparticle on the Fermi 
surface and a,/3 are the spin labels of the paired 
quasiparticles. Fermion statistics requires that the 
pair amplitude obey the anti-symmetry condition, 
A Q/3 (p/) = -Ap a (-p f ). 

Essentially all of the candidates for uncon- 
ventional superconductivity, including the heavy- 
fcrmion and cuprate superconductors, have inversion 
symmetry. This has an important consequence; the 
pairing interaction that drives the superconducting 
transition decomposes into even- and odd-parity sec- 
tors |g), [7). Thus, A a p(pf) necessarily has even or 
odd parity unless there is a second superconducting 
instability into a state with different parity. 

The d-wave model for the cuprates 

Superconductivity in the high T c cuprates is gen- 
erally believed to result from pairing correlations 
within the CuO planes. Knight shift measurements 
below T c indicate that the pairs form spin singlets 
||, and therefore even-parity orbital states. Unlike 
the superfluid phases of 3 He the cuprates exhibit a 
single superconducting phase. The absence of multi- 
ple superconducting phases in the cuprates suggests 
that the orbital pairing state belongs to one of the 
four one-dimensional representations illustrated in 
Fig. 1, for tetragonal crystal symmetry. 
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Fig. 1. Graphical representation the 1-dimensional 
basis functions for even-parity tetragonal supercon- 
ductors. The broken reflection planes and the signs 
of the pairing state are indicated. 

The popular d x 2_ y 2 model for the cuprate su- 
perconductors Q belongs to the B\ g representa- 
tion. This state breaks the reflection symmetries for 
the (110) planes. As a result the order parameter 
changes sign in momentum space as indicated in Fig. 
1, and the excitation gap |A(p/)| has nodes in the 
(110) direction on the Fermi surface. Both features 
lead to novel transport properties in the supercon- 
ducting state at very low temperatures |llj . 

Multi- component Models of UPt^ 

Considerable evidence in support of an uncon- 
ventional pairing state in the heavy-fermion ma- 



terials has accumulated from specific heat, upper 
critical field and various transport measurements, 
all of which show anomalous properties compared 
to those of conventional superconductors (see Refs. 
({3, |l^|) for original references). However, some of 
the strongest evidence for unconventional supercon- 
ductivity comes from the multiple superconducting 
phases of UPt 3 |j| which are a strong indication 
of a multi-component pairing amplitude that neces- 
sarily breaks one or more symmetries of the normal 
state. 

There are several important features of the H-T 
phase diagram, (i) There are two zero-field supercon- 
ducting phases with a difference in transition tem- 
peratures of AT C /T C ~ 0.1. (ii) A change in slope of 
the upper critical field (a 'kink' in H^) is observed 
for Hie. (iii) There are three flux phases, and 
the phase transition lines separating the flux phases 
appear to meet at a tetracritical point for all orien- 
tations of H relative to c. Three basic models that 
have been proposed to explain the phases of UPt3. 
1. Multi- component Order Parameter coupled to a 
Symmetry- Breaking Field (SBF). These are models 
based on a primary order parameter belonging to 
a higher dimensional representation of the symme- 
try group of the normal state. For UPt3, which has 
hexagonal (Dq^) symmetry, there are four 2D rep- 
resentations, £i ff (u) a- n d E2g(u) , with basis functions 
that transform like, 3^Ei± ~ Pz(Px±iPy) for the even- 
parity E\ g representation, and 3^_b 2 ± ~ Pz(j>x±iPy) 2 
for the odd-parity E^ u states. The order parameter 
is then a complex two-component vector 77 = (771, 772) 
that transforms according to the relevant 2D repre- 
sentation, and is related to the pairing amplitude by 



A(p f ) = r)+y El:2 + + r]-y Bl , 2 - 



(2) 



with r/± = rji±ir)2- Multiple superconducting phases 
correspond to different stationary solutions of the 
free energy functional for 77. 

The small splitting of the double transition in 
UPt 3 (AT C /T C ~ 0.1) suggests the presence of a 
small symmetry breaking energy scale and an asso- 
ciated lifting of a degeneracy of the possible super- 
conducting states belonging to the 2D representation 
]l5| , |l6| , [TtJ . The second zero-field transition just be- 
low T c in UPt3, as well as the anomalies observed 
in the upper and lower critical fields, have been ex- 
plained in terms of a weak symmetry breaking field 
(SBF) that lowers the crystal symmetry from hexag- 
onal to orthorhombic, and consequently reduces the 
2D E2 (or Ei) representation to two ID representa- 
tions with slightly different transition temperatures. 



The key point is that right at T c all states belonging 
to the 2D representation are degenerate, thus any 
SBF that couples to ff in second-order, and prefers 
a particular state, will dominate very near T c . At 
lower temperatures the SBF energy scale, AT C , is 
a small perturbation compared to the fourth-order 
terms in the fully developed superconducting state 
and one recovers the results of the GL theory for the 
2D representation with small perturbations to the 
order parameter. 

The phase diagram determined by ultrasound ve- 
locity measurements indicates that the phase bound- 
ary lines meet at a tetracritical point for both H||c 
and H _L c. This has been argued to contradict the 
GL theory based on a 2D order parameter |l8|, 

. The difficulty arises from gradient terms of the 
form, [(D x rn)(D y n 2 )* + (D x ri 2 )(D y rii)* + c.c], that 
couple the two components of the order parameter. 
These terms lead to 'level repulsion' effects in the lin- 
earized GL differential equations which prevent the 
crossing of two H C 2(T) curves, corresponding to dif- 
ferent superconducting phases. This feature of the 
2D model has spawned alternative models, designed 
specifically to eliminate the 'level repulsion' effect 
[ [l9| 53|, and to a more specific version of the 2D 
model coupled to a SBF based on the E 2u represen- 
tation 

2. Accidental degeneracy of two order parameters. 
These models are based on two primary order param- 
eters that are unrelated by symmetry, i.e. belong to 
different irreducible representations of the symmetry 
group, and are accidentally nearly degenerate. The 
model of two ID representations, A 2 and B\ ("AB 
model") is a specific example ^o). The motivation 
behind this model is that by choosing the two repre- 
sentations appropriately one can guarantee that the 
'level repulsion' terms in the Ginzburg-Landau equa- 
tions for H c2 are absent by symmetry |l8[ |2(| . What 
is required is that the two order parameters corre- 
sponding to the two irreducible representations, have 
different signatures under reflection, or parity. In 
this case the second-order gradient coupling between 
the two order parameters vanishes, and the apparent 
tetracritical point is present for all field orientations. 
The drawback is that the accidental degeneracy mod- 
els provide no explanation for the near degeneracy 
of T c , the observed correlation between AT C and the 
AFM order parameter, or the pressure-temperature 
phase diagram. 

3. Enlarged Symmetry Group Models. These are hy- 
brid models that assume an accidental degeneracy 
in the form of a larger symmetry group for the nor- 



mal phase than one expects based on the atomic and 
electronic structure of the crystal. A SBF is then in- 
voked to lift the degeneracy of a higher dimensional 
representation for the larger symmetry group. By as- 
suming a larger group than Dq^ one can again elim- 
inate the 'level repulsion' terms exactly in the GL 
limit by judicious choice of the primary representa- 
tion. Two different versions of the enlarged symme- 
try group model have been proposed; one based on 
an enlarged orbital symmetry group and the other 
based on the assumption of no spin-orbit coupling. 
The orbital model starts from the full rotation group, 
Gspace — 50(3) and produces multiple phases by 
crystal field splitting of the pairing states |^2| . This 
is a flexible model, but there is no strong evidence to 
support treating UPt3 as an isotropic material, even 
approximately; and there is no explanation for the 
correlation between the AFM order and the multiple 
superconducting phases. 

The spin-channel version of this model |l9| as- 
sumes a symmetry group composed of independent 
orbital and spin rotations with Q sp i n — SU(2). The 
multi-component order parameter in this model cor- 
responds to the three spin-triplet amplitudes defined 
in terms of a 3D complex d vector. The 'level repul- 
sion' terms are absent for a ID orbital representa- 
tion. Coupling of the AFM order parameter to the d 
vector is proposed to split the transitions for the dif- 
ferent triplet sub-states. This model requires weak 
spin-orbit coupling, which is at odds with theoreti- 
cal estimates for the spin-orbit coupling energy [fl2| . 
Furthermore, weak spin-orbit coupling is in contra- 
diction with the observation of anisotropic Pauli lim- 
iting p3| ; a spin-triplet order parameter with weak 
spin-orbit coupling would not exhibit Pauli limiting 
for any field orientation (see also Ref. pll ) . 

More on the SBF model for UPt% 

In the absence of accidental near degeneracy, a 
SBF is essential for lifting the degeneracy of the 
pairing states near T c and producing multiple su- 
perconducting phases jl6| [l?], A natural can- 
didate for a SBF in UPt 3 is the AFM order in the 
basal plane |2q] . In this case the GL functional in- 
cludes a coupling of the AFM order parameter to 
the superconducting order parameter; ^-"sbf[??] = 
e M s 2 J d 3 x (hi | 2 - \r} 2 1 2 ) , where M s is the AFM or- 
der parameter and the coupling parameter eM 2 de- 
termines the magnitude of the splitting of the su- 
perconducting transition. The analysis of this GL 
theory, including the SBF, is given in Ref. (]T(i|); 
some of the main results summarized below. 



A double transition is predicted with a splitting of 
T c proportional to AT C oc Af s 2 . Support for the SBF 
model of the double transition comes from pressure 
studies of the superconducting and AFM phase tran- 
sitions. Heat capacity measurements by Trappmann, 
et al. [^6| show that both zero-field transitions are 
suppressed under hydrostatic pressure, and that the 
double transition disappears at ~ 4 kbar. Neu- 
tron scattering experiments reported by Hayden, et 
al. |2^] show that AFM order disappears on the same 
pressure scale, at p c ~ 3.2 kbar. 

The low temperature phase (T < T c „) is pre- 
dicted to have broken T symmetry. This phase is 
doubly degenerate: ff± ~ (a(T) , ±i b(T)) , reflecting 
the two orientations of the internal orbital momen- 
tum of the ground state. 
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Fig. 2. Calculations of and vs - T for 
an odd-parity, triplet pairing state with the d vector 
locked to the c direction ]28j . The inset shows the ra- 
tio H^/H^2 vs. T/T c . Note that h]} 2 is suppressed 
by paramagnetism (/x ~ /j,b) at low temperatures, 
while H^2 is independent of the paramagnetic cou- 
pling. The data of Shivaram, et al. [^3[ are shown 
as the open and closed circles. 

Another key feature of the UPt3 phase diagram is 
the unusual anisotropy of H C 2 shown in Fig. 2. The 
low-temperature anisotropy is explained in terms of 
anisotropic Pauli limiting of an odd-parity, spin- 
triplet state with the d-vector parallel to the c di- 
rection pg|| . This restricts one to the Yj2 U or 
symmetry classes among the four possible 2D repre- 
sentations. 

The SBF is essential for producing an apparent 
tetracritical point, and at a semi-quantitative level, 
can account for the magnitudes of the slopes near the 
tetracritical point pi], p9|. However, there are open 



questions regarding both the nature of the AFM or- 
der and its coupling to superconductivity fl30| |. As 
this model for UPt3 illustrates, GL theory provides 
a central link between experiments, in this case the 
phase diagram, and the more microscopic Fermi- 
liquid theory which we discuss below. 

2. Fermi-Liquid Theory of Heavy Fermion 
Superconductivity 

Conduction electrons in metals interact strongly 
with each other and with the lattice. These interac- 
tions lead to correlations among the electrons, and 
we have to view conduction electrons, in general, as a 
system of correlated fcrmions. A particularly impor- 
tant system of strongly correlated electrons are the 
conduction electrons of heavy fermion metals. There 
is strong evidence that these heavy electrons form a 
Landau Fermi liquid at low temperatures, and that 
their superconducting states are well described by 
the Fermi liquid theory of superconductivity. In the 
following we give a brief introduction to this the- 
ory, then discuss recent applications to transport in 
heavy fermion superconductors with unconventional 
pairing. 

Landau showed that an ensemble of strongly in- 
teracting fermions may be described by a distribu- 
tion function for quasiparticle excitations, and that 
this distribution function obeys a classical transport 
equation, the Boltzmann-Landau transport equa- 
tion. It took more than 10 years after Landau's the- 
ory of normal Fermi liquids, and the breakthrough 
in the theory of superconductivity by BCS, to es- 
tablish a complete Fermi-liquid theory of supercon- 
ductivity. Earlier general theories of superconduc- 
tivity lacked the quasiclassical aspects of Landau's 
Fermi-liquid theory. The first complete quasiclas- 
sical (QC) theory of superconductivity was formu- 
lated in a series of publications by Eilenberger pl[ , 
Larkin and Ovchinnikov |3^| and Eliashberg [[53 . It 
is presented and discussed in several review articles 
HJ, HI The quasiclassical theory allows one to 
calculate all superconducting phenomena of interest, 
including transition temperatures, excitation spec- 
tra, Josephson effects, vortex structures, the electro- 
magnetic response, etc. In this theory the dynami- 
cal degrees of freedom of electronic quasiparticles are 
described partly by classical mechanics, and partly 
by quantum statistics. The classical degrees of free- 
dom are the motion in p-R phase space; i.e. quasi- 
particles move along classical trajectories. Quantum 
degrees of freedom are the spin of a quasiparticle 
and the particle-hole degree of freedom, which form 



a four-dimensional Hilbert space of "internal degrees 
of freedom" . 
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Fig. 3. Leading order self-energy diagrams of 
Fermi-liquid theory. The vertices (shaded circles) 
represent the sum of all high-energy processes and 
give rise to interactions between the quasiparticles 
(smooth propagator lines), phonons (wiggly propa- 
gator lines) and impurities (dashed lines) . The order 
in the parameter 'small' ~ e/Ef is indicated for each 
diagram. 

Derivations of the Boltzmann-Landau transport 
equation from first principles |37| |3£| use many- 
body Green's function techniques, and lead to ex- 
plicit expressions for the various terms of the trans- 
port equation in terms of self-energies. The self- 
energies describe the effects of electron-electron, 
electron-phonon and electron-impurity scattering. 
The complete set of diagrams for the leading or- 
der self-energies are shown in Fig. 3; the filled 
circles are vertices representing the effective lat- 
tice potential, which determines the quasiparti- 
cle Fermi surface and Fermi velocities, as well 
as quasiparticle-quasiparticle, quasiparticle-phonon, 
and quasiparticle-impurity interactions. We follow 
Landau and consider these vertices as phenomeno- 
logical parameters of the Fermi-liquid model, but in 
principle they can be obtained from the full many- 
body theory. The selection of leading order self- 
energies holds for the normal and the superconduct- 
ing states. The vertices are not affected in leading 
order by the superconducting transition; supercon- 
ductivity affects only the electron propagators. The 
special achievement of Eilenberger, Larkin, Ovchin- 
nikov, and Eliashberg was to convert Dyson's equa- 
tions for the electron propagators into transport 
equations for quasiclassical propagators. The result- 
ing equations define the quasiclassical theory of su- 
perconductivity. 



The central equation of the quasiclassical the- 
ory is the transport equation for the quasiclassical 
Keldysh propagator, g K , 

(ef 3 - v- a R ) <g> g K - g K <g> (ef 3 - v - a A ) - 

ihv r Vg K = 0. (3) 



r K 



• 9 + 9 



The retarded, advanced, and Keldysh propaga- 
tors, g R ' A ' K (pj t R; e, t), as well as the self-energies, 
a R,A ' K (pf, R; e, i), and the external potentials, 
u(p/,R, t), are 4x4 Nambu matrices, acting on the 
4D Hilbert space of internal degrees of freedom. The 
<g>-product stands for the usual 4 x 4-matrix prod- 
uct and a product in the energy-time variables. De- 
tails of this compact notation are explained in the 
reviews of the quasiclassical theory |34|, |3(J . The 
determination of g K from the transport equation re- 
quires knowledge of the external potentials, the ad- 
vanced, retarded and Keldysh self-energies, and the 
advanced and retarded quasiclassical propagators. 
These propagators are auxiliary quantities which in 
general have no direct physical interpretation, except 
in the adiabatic limit |34|j where they determine the 
local quasiparticle density of states. The retarded 
and advanced propagators are solutions of 

[ e r 3 - e - a R - A , g R ' A ] g + ihv, ■ Vg R > A = . (4) 

The physically relevant set of solutions of Eqs. (||) 
and (0) must satisfy the normalization conditions, 



g R > A ®g R > A = 



Vi, g R ®g K + g K ®g A 



0. (5) 



Measurable quantities such as the charge current 
density, j (R, t) , can be calculated from the diagonal 
components of the quasiclassical propagator, g . For 
example, the charge current density is given by 

j(R> * } = / jSff^Tl I ^ ev/ gK(ph R; e ' * } ' 

(6) 

The quasiclassical equations (||)-(||) are supple- 
mented by self-consistency equations for the the qua- 
siclassical self-energies. These equations are shown 
in diagrammatic notation in Fig. 3. They include, 
for example, the "gap equation" , which is the self- 
consistency equation for the off-diagonal self-energy. 
Explicit forms of the quasiclassical self-energies can 
be found in Refs. g, [3§, §§]. 

The quasiclassical theory is especially well suited 
for studying unconventional superconductors. Some 
of the results of the theory for conventional s-wave 
superconductors are simply generalized to unconven- 
tional superconductors. For example, the density 



of states N(pf, e) of a homogeneous spin-singlet su- 
perconductor in equilibrium has the standard BCS 
form, N(pf,e) = N f $l[e/y/\ A | 2 -e 2 ], but with the 
isotropic gap |A| replaced by an anisotropic gap, 
|A(p/)|, for each point on the Fermi surface. 

However, inhomogeneous and non-equilibrium 
situations exhibit more striking differences between 
conventional and unconventional superconductors 
that reflect both the broken symmetries of the or- 
der parameter and the coherence properties of the 
superconducting state. 

Particle-Hole Coherence in QC Theory 

Of special importance for understanding super- 
conducting phenomena are the quantum-mechanical 
"internal degrees of freedom", the spin and the 
particle-hole degrees of freedom of an electron. 
Quantum coherence between particle excitations and 
hole excitations is a key feature of the BCS theory of 
superconductivity and the origin of all non-classical 
effects in superconductors (e.g. supercurrents, co- 
herence factors in transition amplitudes, Andreev 
reflection, etc.). Particle- hole coherence is incorpo- 
rated into the quasiclassical theory by grouping par- 
ticle excitations (occupied one-electron states above 
the Fermi energy) and hole excitations (empty one- 
electron states below the Fermi energy) into a dou- 
blet. For clean superconductors, and for long wave- 
length spatial variations, £o = hvf/2nT c ^> h/pf, 
we can make a quasiclassical envelope approximation 
for the Bogoliubov amplitudes; the resulting particle 
and hole amplitudes obey Andreev's equation, [[f2| 

(ef 3 - A( P/ , R)) PJ + ihv f • V <p Pf = , (7) 

where tp Pf = {u Pfl v Pf ) are the quasiclassical par- 
ticle (u) and hole (v) amplitudes. Andreev's equa- 
tion is a first-order differential equation for excita- 
tions propagating along straight-line trajectories de- 
termined by the Fermi velocity v/(p/) at a point p/ 
on the Fermi surface. In the superconducting state, 
the order parameter A( P /,R) mixes the normal- 
state excitations coherently into two branches of 
particle-like (p/ • v + > 0) and hole-like (p/ • v_ < 0) 
excitations. 

Coherence between particle and hole excitations 
leads to dramatic effects on the excitation spectrum 
of an unconventional superconductor in the vicin- 
ity of an impurity or a surface. Consider an exci- 
tation that is incident on a specular surface from 
the bulk region of a superconductor with broken re- 
flection symmetry perpendicular to the plane of the 



interface. When an excitation reflects off the surface 
elastically its momentum shifts to a new point on 
the Fermi surface, p/ — » p^.. Thus, the incident and 
reflected wavepacket propagate through different or- 
der parameter fields, A(p/,R) vs. A(p^,R). As a 
result, surface scattering generally leads to Andreev 
scattering, a process of "retro-reflection" in which a 
particle-like excitation undergoes branch conversion 
into a hole-like excitation with reversed group veloc- 
ity. Bound states may occur at energies for which 
the phases of multiply-reflected particle- and hole- 
like excitations interfere constructively. 

The effects of a surface on particle-hole coher- 
ence is most pronounced if the scattering induces 
a change in sign of the order parameter along the 
classical trajectory. This occurs for a d x 2_ y 2 super- 
conductor with a (110) surface, an E lg or E 2u order 
parameter and a surface normal to the c axis. If 
a sign change of the order parameter occurs along 
the trajectory, then a zero-energy bound state forms 
at the surface with equal amplitudes for the particle 
and hole components ]3j| . These states are expected 
to give rise to zero-bias anomalies in the conductance 
for NIS tunnel junctions (4fj, |4l], [TJ) . 

Impurity-induced Andreev bound states 

The novel effects of impurities in unconventional 
superconductors, including pairbreaking, can be un- 
derstood in similar terms. Fermion bound states 
are formed by impurity-induced Andreev scattering. 
Figure 4 illustrates the connection between potential 
scattering and the development of an Andreev bound 
state for s-wave impurities in the d x i _ y 2 model; sim- 
ilar physical processes apply to unconventional mod- 
els of the heavy fermion superconductors. Impurities 
give rise to elastic scattering of states near the Fermi 
surface, with transition matrix elements given by 
u(p/,pj). The scattering amplitude is determined 
by the t-matrix, 

F' A ( P/ , p/, R; e, t) = u(p f ,p/)+N f J dk f 

u(p f , k f )g R > A (k f ,R; e, t) ® t R < A (k f , p/, R; e, t) (8) 

The i-matrix sums two types of repeated scattering 
processes to all orders: (i) potential scattering pro- 
cesses with a change in momentum, but no change in 
the internal state, and (ii) Andreev scattering pro- 
cesses, i.e. branch conversion with no change in mo- 
mentum, but reversal of the group velocity. The An- 
dreev processes are induced by the sign changes of the 
order parameter that result from scattering around 



the Fermi surface as shown in Fig. 4. Bound states 
also form from impurity-induced Andreev scattering, 
with an energy that depends on the scattering phase 
shift. Neglecting broadening effects from the contin- 
uum states above the gap, the bound state energy 
is given by thound = Ao cos^o. However, the contin- 
uum excitations are nearly gapless in the vicinity of 
the nodes of A(p/), so bound states at finite energy 
are broadened into resonances. For strong scattering 
the bound state occurs near zero energy and the res- 
onance width is narrow; e.g. in the unitarity limit 
the impurity bound states appear at e = and are 
sharp. 




U(P,P r ') =u„ 




Fig. 4. Scattering of quasiparticles by an impu- 
rity induces Andreev scattering, and for unconven- 
tional pairing states the formation of Andreev bound 
states or resonances. 
In a dilute alloy a finite density of impurities leads 
to a finite density of impurity states, an Andreev 
band The impurity bandwidth and density of 

states can be calculated from the self-consistent t- 
matrix and the leading order impurity self-energy 
terms (diagrams 3e). The bandwidth is given by 

(7(|A(P/)| 2 +7 2 r 1/2 ) 



cot 2 <5 + (7(|A( P/ )| 2 +7 2 ) 



where r„ = n^p/nNf is the maximum impurity 
scattering rate in the the normal state for a fixed 
concentration, ni mp , of impurities and (...) is an av- 
erage over the Fermi surface. In the unitarity limit 
the impurity bandwidth is 7 ~ yj hY u /2Aq ~ ^/7ij mp . 

For 7 < e <C Ao the low energy spectrum is 
dominated by continuum excitations in the vicinity 
of the nodal lines; for the E\ g and E 2u models of 
UPt 3 N(e) w Nf(e/A ). Below e w 7 the impurity 
band dominates with N(e) ~ N(0) s» N f (-f/A ). 
These two contributions to the spectrum also give 
rise to different features in the transport proper- 
ties of the heavy fermion superconductors (l(], ^| 
To calculate the transport coefficients in the super- 
conducting state we need the Keldysh propagator for 
these branches of the excitation spectrum. Linear re- 
sponse functions for electrical and thermal transport 
in strong-coupling, dirty and unconventional super- 
conductors are derived from the quasiclassical trans- 
port equations in Refs. [prjl llOl pi. 



3. Thermal Conductivity of UPt 3 - E 2u Model 

The anisotropic thermal conductivity coefficients 
for the Ei u model are shown in Fig. 5 as a function of 
temperature, and compared with the data of Lussier, 
et al. Q for T > 0.1T C . The fits of the E 2u model 
to the low temperature region, 0.1 < T/T c < 0.5 
are very good for both directions of heat flow. More 
detailed analysis shows that both the E 2u and E\ g 
pairing states give excellent fits to the thermal con- 
ductivity for T > 0.1 T c , and that ultra- low tem- 
perature (i.e. T < 7) heat transport measurements 
and impurity studies should distinguish these pairing 
states [ fj"i"l . 



1.4 




1.2 



Fig. 5. The normalized thermal conductivity of 
the E2 U ground state compared with the data 
of Ref. Ma|. The impurity scattering rate is To = 
O.Ol7rT c o and the phase shift is 80 = 90°. The slope 
and curvature parameters are fii — 2.0 and fi2 = 4.0. 

Ultra-Low Temperature Heat Transport 

Signatures the pairing symmetry arise from both 
the low-energy continuum states and the impurity- 
induced Andreev band. This novel metallic band 
deep in the superconducting state gives rise to uni- 
versal transport coefficients for fc^T <C 7 JlO| . 
Whether or not a universal limit develops at low tem- 
peratures depends sensitively on the nodal structure 
of the excitation gap. Both the Ei s and F,2 U models 
for the order parameter of UPt3 lead to an excitation 
gap with a nodal line in the basal plane and point 
nodes along the ±c directions jl7|. The difference in 
the excitation spectrum for these two states is that 
the gap opens linearly for small angles away from the 
c direction for the Ei fl order parameter, but quadrat- 
ically for E 2u order parameter; |Ae 1js (i?)| — Ao/ii|i?| 
and |A Eaa (0)| ^ A ^ 2 |^| 2 - This difference is re- 
flected in limT^o k c /T, which is universal for the E2 U 



model, but non-universal for the Ei 9 model |11 



K c /T ' 
K c /T c 



—N 2 v 2 f 
3 / / 



-Njv 



f u f 




(10) 



(11) 



For heat flow in the basal plane the thermal con- 
ductivity is determined by the line node, and is uni- 
versal for both E\ g and Em pairing states. The lead- 
ing temperature corrections to k/T are given by a 
Sommerfeld expansion for T < 7 < T c . The T 3 cor- 
rections are non- universal and in the unitarity limit, 
7 2 oc rii m p, so the coefficient of the T 3 term scales 
as f/ni m p for the E2„ state and 1 / y/n- lmp for the Ei g 
state. Observed universality for T <C 7 and the scal- 
ing of the T 3 correction with 7ii mp would be strong 
tests of the symmetry of the order parameter. 

We thank M. J. Graf, D. W. Hess, and S.-K. Yip 
for many contributions to the work reported here. 
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